Antiferromagnetic domain walls in lightly doped layered cuprates 
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Recent ESR data shows rotation of the antiferromagnetic (AF) easy axis in lightly doped layered 
cuprates upon lowering the temperature. We account for the ESR data and show that it has 
significant implications on spin and charge ordering according to the following scenario: In the high 
temperature phase AF domain walls coincide with (110) twin boundaries of an orthorhombic phase. 
A magnetic field leads to annihilation of neighboring domain walls resulting in antiphase boundaries. 
The latter are spin carriers, form ferromagnetic lines and may become charged in the doped system. 
However, hole ordering at low temperatures favors the (100) orientation, inducing a 7r/4 rotation in 
the AF easy axis. The latter phase has twin boundaries and AF domain walls in (100) planes. 

PACS numbers: 75.25.+z,75.60.Ch,75.80.+q,61.72.Mm 



The effect of holes on antiferromagnetism is a cen- 
tral issue in high temperature superconductivity. In 
a remarkable recent experiment the role of doping in 
Y\- x Ca x Ba,2CuzO§ [x ~ 0.008) was studied by Janossy 
et ali using an ESR technique. The data shows that the 
antiferromagnetic (AF) polarization rotates from [100] at 
high temperatures to [110] around » 40K. A correspond- 
ing anomaly was also seen by /iSR in Y\_ x Ca x Ba,2Cu?,0§ 
as well as in La,2- x Sr x CuOi compounds*. In the un- 
doped compound (x = 0) the polarization remains in the 
[100] directions down to low temperatures. 

The ESR data* also shows that AF domains are 
present and that their relative intensity is controlled by 
magnetic fields. In the high temperature phase with field 
in the [100] direction, the [010] polarized domains are 
preferred since the canting of AF spins has a higher sus- 
ceptibility. This defines a " depinning" field at which the 
unpreferred domains are diminished. E.g. at high tem- 
peratures and above a field of w IT 80% of the domains 
are polarized perpendicular to the field. The ESR data 
shows, curiously, that the depinning field of the [110] AF 
of Yi- x Ca x Ba2CuzOQ (with field in the [110] direction) 
is higher than that of the [100] AF in YBa 2 Cu 3 6 (field 
in the [010] direction), both at low temperatures. 

In the present work we suggest that the transition in 
the AF polarization correlates with a lattice distortion 
and with condensation of the holes into a charge density 
wave (CDW). Our reasoning involves a scenario for spin 
and charge ordering with the following steps: (i) At high 
temperatures and weak magnetic fields the [100] polar- 
ized AF leads to an orthorhombic structure; the AF do- 
main walls (DW's) coincide then with twin boundaries at 
(110) planes, (ii) Strong magnetic fields favor one of the 
AF domains so that DW's annihilate. We show that this 
process results in antiphase boundaries (APB's) which 
consist of ferromagnetic lines, (iii) At low temperatures 
hole ordering can occur by binding states localized at 
the APB's. However, hole ordering favors the (100) ori- 
entation, i.e. a CDW with maxima on (100) planes and 
an (incommensurate) wavevector in the [100] direction. 
This CDW favors an AF polarization in a diagonal [110] 
direction which is weakly orthorhombic with (100) twin 



boundaries and DW's. 

An AF polarized in the [100] axis favors an orthorhom- 
bic phase since dipoles interact differently for — > <— and 
1 1 pairs. This magnetoelastic coupling is expected to be 
rather small in the cuprates, leading to a change in lattice 
constant of order— 10~ 5 — 10~ 6 . In principle ESR analysis 
can detect the resulting change in crystal field*, however, 
the effect seems too weak in the case AF Y Bcl2CuzOq. 
Yet, this coupling can be sufficiently strong in samples 
of size w 1mm so that twin boundaries are necessary to 
relieve macroscopic strains. 

We proceed to present evidence for steps (i) , (iii) and 
then derive step (ii). We claim that an orthorhombic 
structure is in fact supported by the mere observation 
of AF domains. The presence of AF domains and the 
associated domain walls* are well known to result from 
either entropy effects (when temperature is fairly close 
to the Neel temperature), or from disorder or from a 
magnetoelastic coupling. The ESR experiment* has used 
high purity samples and temperature was well below the 
Neel temperature. Hence the presence of DW's is a strong 
indication for the presence of the magnetoelasic coupling, 
supporting step (i) of the scenario above. 

Step (iii) is consistent with data on the compounds 6 
Lai.s75Bao.i25-xSr x Cu04,. These compounds maintain 
fixed carrier density while allowing for both tetragonal 
and orthorhombic phases. In the tetragonal phase an in- 
commensurate CDW (as well as an incommensurate spin 
density wave) with wavevector in the [100] direction is 
present. For systems near the orthorhombic boundary 
(less orthorhombic phase) the CDW wavevector deviates 
slightly from the [100] axis and finally the CDW disap- 
pears in the orthorhombic phase. 

We propose then that the orthorhombic structure in- 
duced by the [100] polarized AF structure is unstable 
when holes are added and form a [100] CDW. The AF 
polarization then rotates to the [110] direction which, as 
shown below, is weakly orthorhombic and has (100) twin 
boundaries. We proceed now to analyze an effective free 
energy and study various domain walls. In particular step 
(ii) is shown, i.e. the collapse of two neighboring mag- 
netic DW's (e.g. by applying a magnetic field) leads to 



2 




FIG. 1: Twin boundary with spin polarizations (arrows) ex- 
hibiting an AF domain wall. The dashed line is in a (110) 
plane and is the s = line defining a TB position. The dis- 
placement u(s) is parallel to this line and yields the strain 
e = du/ {y/2ds) interpolating between the two orthorhombic 
variants. 



antiphase boundaries which can be detected by a variety 
of experiments. 

Consider first the elasticity part, i.e. the tetragonal- 
orthorhombic (T-O) transition and formation of twin 
boundaries (TB's). The orthorhombic distortion, for 
weak strain, corresponds to a tetragonal lattice displaced 
in the [110] direction with a displacement u that is linear 
in the coordinate s in the [110] axis, i.e. u ~ ±s (see 
Fig. 1). The ± signs correspond to the two variants of 
the orthorhombic phase. The free energy, in terms of the 
strain e — du/(\/2ds), has then the form 7 

where Ti{e) has a double minima corresponding to the 
two variants. A static solution interpolating smoothly 
between the two variants is then possible 7 *, and is de- 
scribed by a function €tb{s)- This function has neces- 
sarily an s value where £tb( s ) — which defines a twin 
boundary location, hence the TB is on a (110) plane. 
In general, a TB orientation is uniquely determined by 
the type of structural transformation, i.e. smooth inter- 
polation between variants across other planes leads to 
diverging elastic energies^. The formulation of Eq. 
is a simple demonstration that for the T-0 transition 
TB's are on (110) planes. A macroscopic strain which 
maintains the overall sample shape, e.g. by presence of 
a parent (tetragonal) phase 9 , imposes a TB array and 
determines its periodicity. 



We proceed now to study the T-0 transformation in- 
duced by a [100] polarized AF as well as the possible 
AF domain walls coexisting with TB's. An AF domain 
wall is defined by a localized rotation of the polariza- 
tion angle 9 (e.g. relative to the [100] axis) by ±7r/2, 
i.e. neighboring domains are polarized along [100] and 
[010], respectively. The coupled AF and strain structure 
is illustrated in Fig. 1. Note that antiparallel spins favor 
shorter bonds, therefore a TB induces a ir/2 rotation in 
the AF polarization, as shown in the figure. 

We claim that neighboring DW's have the same rota- 
tion angle ir/2 (or —ir/2), rather than opposite 7r/2 and 
— 7t/2. To show that the latter has no solution note first 
that the macroscopic strain which imposes the TB array 
does not affect the local differential equation for 9(s), 
instead it determines the average TB spacing 9 being an 
integration constant. Hence the existence of a static so- 
lution 9(s) can be deduced from dynamical stability con- 
siderations of the AF structure by itself. In particular 
DW's with opposite 7r/2 and — 7r/2 rotations are unsta- 
ble since by approaching each other they can annihilate 
and form the lower energy ground state; in contrast, two 
7r/2 DW's when approaching each other form a higher 
energy antiphase boundary (see below), hence they are 
stable and a static solution is possible. 

To illustrate this idea more explicitly we consider a 
slowly varying magnetization with a magnetoelastic cou- 
pling energy of the form — ae cos(2#); in particular a [110] 
polarization (9 — 7r/4) does not couple to this strain (the 
7r/4 state has a weaker coupling to a different strain, as 
considered below). F\ in Eq. J2J is replaced now by a 
quadratic term + \be 2 since the T-0 transition is driven 
by the magnetoelastic coupling, hence the full free energy 
is 

« « + ^ 3+ H!) 2+ HS) 2 (2) 

where the constants b, c, d are positive and the last term 
represents the spin stiffness. (A crystal field ~ cos(4#) 
from spin orbit coupling is neglected; in fact it must be 
small to ensure that 9 — is a ground state.) The two 
orthorhombic variants correspond to 9 = 0, 7r/2 with the 
strain e = ±a/b. 

Domain walls are solutions of the minimum condition 
for Eq. 0, 

<9 2 e 

- acos(29) + be - c—r- = 
o z s 

d 2 9 

-2aesin(26>) + d—r- =0. (3) 

o z s 

A TB which interpolates between the variants e = ±a/b 
necessarily leads to rotation of the AF polarization by 
7r/2, as in Fig. 1. We show next that a periodic TB 
array for which 9(s) is a monotonic function is a valid 
solution. Consider a segment of such a TB array where 
9 = at s = —so, becomes 9 = ir/4 at s — and finally 
rotates to 9 = n/2 at s = sq; in the same interval the 
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FIG. 2: The 1st line shows three AF variants along [100], 
separated by two DW's (with the diagonal polarizations) . The 
2nd line shows the result of annihilation of the two DW's - the 
spins marked by ? continue either left or right domains and 
an antiphase boundary results somewhere in the ? marked 
range. 
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FIG. 3: An APB along a (11) line, centered on the dots, 
showing a ferromagnetic line. Note that the ferromagnetic 
polarization is parallel to that of the AF. 



strain is antisymmetric, interpolating from e near a/b > 
(assuming a > 0) to e near —a/b < 0. Integration of Eq. 
© yields 

d6(s) d9(s). 2a f s , , . . „, S1 , 

In the interval — sq < s < sq sin(2(9(s)) > is symmet- 
ric while e is antisymmetric starting from positive values, 
hence the right hand side of Eq. (0} is positive and ap- 
proaches zero at sq. At the next TB, e is antisymmetric 
starting from negative values while sin(2#(s)) is symmet- 
ric and negative, hence again the right hand side of Eq. 
(@J is positive. We conclude that d9/ds > is a con- 
sistent solution. Furthermore, in the limit c <C d the 
gradients are dominated by the spin stiffness while the 
strain follows the local AF angle, e w (a/6)cos(2#) and 
Eq. © becomes -(a 2 /b) sin(40) + dd 2 9/d 2 s = 0. The 
latter is the sine-Gordon equation with well known peri- 
odic domain structure satisfying 86/ ds > 0. 

A 7r/2 DW is therefore followed by the same sign +7r/2 
DW. This result is, in fact, a robust consequence of the 
energy argument above requiring DW's stability, hence 
it is valid even with additional terms in Eq. ©. A re- 
markable consequence of this analysis is that when two 
neighboring TB's annihilate the result is not a uniform 
ground state, but rather a change of 9(s) by 7r which is 
an antiphase boundary (APB), i.e. all spins on one side 
of the boundary are reversed relative to the ground state 
orientations. The annihilation process of DW's is shown 
in Fig. 2. The 1st line shows the spin arrangement of 
three variants separated by two DW's with the polariza- 
tion angle increasing monotonically. In the 2nd line the 
DW's are eliminated, i.e. the spins are rotated so as to 
be aligned with either the left or the right domain. The 
diagonal spins join their nearest domain, but the hori- 
zontal ones can choose either domain (indicated by a ? 
mark), with any choice resulting in an APB somewhere 
in the ? marked region. Note that the APB carries an 
additional spin 1/2 relative to the AF ground state. 

APB's involve a major rearrangement of the electronic 
coordinates and can result in a localized state, as shown 
by a mean field study 10 , analogous to studies on soli- 
tons in polyacetylene. An APB involves therefore atomic 
scales and its width is of that order. In contrast, an AF 



DW as well as a TB, have a much larger width and there- 
fore lower formation energies. The high energy cost of an 
APB, relative to that of two DW's, provides the stability 
of the DW array with same sign DW's. 

DW's can be eliminated by a strong magnetic field as 
was demonstrated by the experiment of Janossy et ali, 
leading to a single variant. The process can generate 
APB's assuming that they are sufficiently apart and do 
not annihilate each other. APB's carry spin 1/2 per Cu 
along the boundary, i.e. a ferromagnetic line along [110] 
as shown in Fig. 3. These lines are AF in the [001] di- 
rection, yet, the interlayer coupling is weak so that it is 
relatively easy to flip magnetization by a field parallel 
to the AF polarization. Hence, when the magnetic field 
overcomes the interlayer coupling a nonlinear longitudi- 
nal susceptibility \\\ appears due to aligned ferromag- 
netic lines. Note that by rotating the APB of Fig. 3 into 
a (100) plane its lines (within the layers) are AF rather 
then ferromagnetic. 

Consider next an alternative AF phase with polariza- 
tion along [110], stabilized by terms beyond those in Eq. 
(J2J, e.g. CDW ordering (see below). The coupling of 
nearest neighbors along [100] or [010] is now identical 
and orthorhombicity is not induced, i.e. the first term in 
(J5J) vanished at 9 — ir/4. However, second nearest neigh- 
bors form either | , "f or — >, — > pairs which are nonequiv- 
alent. Hence a tetragonal unit cell which is formed by 
the vectors [110], [110], [001] of the original lattice tends 
to contract or expand along [110] and [110], respectively. 
This unit cell is rotated by 7r/4 relative to the previous 
one, hence the axis s' along which an orthorhombic dis- 
tortion alternates is now in the [100] direction, therefore 
twin boundaries that interpolate between 9 = 7r/4 and 
9 = 3ir/4 are in (100) planes of the original tetragonal 
lattice. The detailed description of these TB's is more 
involved than the previous ones since near 9 « 7r/2 the 
previous strain will couple. However, following the sta- 
bility argument above, we expect again a 7r/2 DW to 
be followed by the same sign +7r/2 DW, hence a strong 
magnetic field in a [110] direction would coalesce these 
TB's and lead to (100) APB's. 

The final ingredient in our model is that the added 
holes form a CDW with wavevector along [100] as clearly 
seen in the Lai ^-j^BaQ.\25-xSr x CuOi compounds 6 . In- 
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TABLE I: The table lists topological structures in the AF ordering in lightly doped Cuprates, as well as pertinent experimental 
signatures at various magnetic fields. High temperature corresponds to the [100] polarized AF phase and low temperature 
corresponds to the [110] polarized AF with a [100] CDW phase. 





high temperature 


low temperature 


strong [100] field 


(110) neutral APB 
ferromagnetic lines, nonlinear \\\ 


(100) DW and (100) charged APB 


strong [110] field 


(110) DW 


(100) charged APB 
higher Tcdw 


no field 


(110) DW 


(100) DW and (100) charged APB 
stronger depinning than for the undoped case 



deed the charged APB in the mean field calculation 10 
is more strongly bound when the APB is along a (100) 
plane. These (100) APB's form a periodic array of stripes 
equivalent to a CDW with wavevector in the [010] direc- 
tion. In a [100] polarized AF these (100) APB's would 
cross (110) DW's or (110) APB's; the [100] CDW there- 
fore favors a rotated [110] polarized AF so that crossings 
with DW's are avoided. Furthermore, in a strong [110] 
field the [110] polarized AF has the appropriate ingredi- 
ents for forming stripes, namely (100) APB's, facilitating 
the CDW formation. We expect therefore that the tran- 
sition temperature Tcdw into a CDW be higher with 
a field in this [110] direction. The predictions for DW's 
and APB's are summarized in table I, as well as pertinent 
experimental signatures. 

We show now that our model accounts also for the ob- 
servation that depinning field in the [110] polarized AF 
Y\- x Ca x Ba<2.Cu-sO§ is larger than that of the [100] po- 
larized AF in YBa2CusOe, both at low temperatures^. 
Depinning involves annihilation of DW in pairs leading 
to a single domain. In Y\- x Ca x Ba2Cuj,0§ some of the 
DW need to cross a charged APB, a process which has 
a barrier. In Y BaiCu-iO§ the process involves TB an- 
nihilation which is continuously achieved as field is in- 
creased with no barrier. Hence a larger depinning field 
in Yi- x Ca x Ba 2 Cu 3 O e . 

Finally, we propose a variety of experiments for ob- 



serving APB's and probing our scenario: (i) APB's are 
expected to have intragap states^, hence optical absorp- 
tion should show new lines at high magnetic fields, (ii) 
The ferromagnetic nature of (110) APB lines (Fig. 3) 
results in a nonlinear longitudinal susceptibility x\\- (iii) 
By adding holes the APB intragap state can be charged, 
leading to a spinless charge carrier. These charges affect 
the optical absorption in (i) as well as the susceptibility in 
(ii). (iv) The CDW onset is facilitated by a strong mag- 
netic field in the [110] direction, hence Tcdw is higher 
in this case. 

In conclusion, our scenario accounts for the unusual 
ESR datai and predicts a variety of topological struc- 
tures (table I). We propose that experiments based on 
these observations can supplement the ESR data as well 
as clarify a central issue in high temperature supercon- 
ductivity - the nature of doping in AF layered cuprates. 
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